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ABSTRACT 

Context. We present the analysis of the baryonic content of 52 X-ray luminous galaxy clusters observed with Chandra in the redshift 
range 0.3 - 1.273. 

Aims. Our study aims at resolving the gas mass fraction in these objects to place constraints on the cosmological parameters Qm, 
and the ratio between the pressure and density of the dark energy, w. 

Methods. We deproject the X-ray surface brightness profiles to recover the gas mass profiles and fit a single thermal component to the 
spectrum extracted from a region around the cluster that maximizes the signal-to-noise ratios in the observation. The measured values 
of the gas temperature are used to evaluate the temperature profile with a given functional form and to estimate the total gravitating 
mass in combination with the gas density profiles. These measured quantities are then used to statistically estimate the gas fraction 
and the fraction of mass in stars. By assuming that galaxy clusters are representative of the cosmic baryon budget, the distribution 
of the cluster baryon fraction in the hottest (Tgas > 4 keV) systems as a function of redshift is used to constrain the cosmological 
parameters. We discuss how our constraints are affected by several systematic effects, namely the isothermality, the assumed baryon 
fraction in stars, the depletion parameter and the sample selection. 

Results. By using only the cluster baryon fraction as a proxy for the cosmological parameters, we obtain that ilm is very well 
constrained at the value of 0.35 with a relative statistical uncertainty of 11% (la level; w — —1) and a further systematic error of 
about (—6, +7)%. On the other hand, constraints on Qa (without the prior of flat geometry) and w (using the prior of flat geometry) 
are definitely weaker due to the presence of greater statistical and systematic uncertainties (of the order of 40 per cent on Qa and 
greater than 50 per cent on w). If the WMAP 5-year best-fit results are assumed to fix the cosmological parameters, we limit the 
contributions expected from non-thermal pressure support and ICM dumpiness to be lower than about 10 per cent, also leaving room 
to accommodate baryons not accounted for either in the X-ray emitting plasma or in stars of the order of 1 8 per cent of the total cluster 
baryon budget. This value is lowered to zero for a no-flat Universe with f^A > 0.7. 

Key words, galaxies: cluster: general - galaxies: fundamental parameters - intergalactic medium - X-ray: galaxies: clusters - 
cosmology: observations - dark matter. 



1. Introduction 

Several tests have been suggested to constrain the geometry and 
the relative amounts of the matter and energy constituents of the 
Universe (see Peebles & Ratra 2003 and references therein). A 
method that is robust and complementary to others is obtained 
using the gas mass fraction, /gas = Afgas/Mtot, as inferred from 
X-ray observations of clusters of galaxies. In this work, we con- 
sider two independent methods for our purpose: (i) we compare 
the relative amount of baryons with respect to the total mass ob- 
served in galaxy clusters to the cosmic baryon fraction to provide 
a direct constraint on fl^ (this method was originally adopted by 
White et al. 1993 to show the limitation of the standard cold dark 
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matter scenario in an Einstein-de Sitter Universe), (ii) we limit 
the parameters that describe the geometry of the universe by as- 
suming that the gas fraction is constant in time, as first suggested 
by Sasaki (1996). 

Starting from these pioneering works, many studies have fol- 
lowed this approach to constrain the cosmological parameters 
(White & Fabian 1995, David et al. 1995, Ettori & Fabian 1998, 
Rines et al. 1999, Roussel et al. 2000, Allen et al. 2002, Ettori, 
Tozzi & Rosati 2003, Castillo-Morales & Schindler 2003, Sadat 
et al. 2005, Allen et al. 2008). We present here a revised and 
updated version of the work discussed in Ettori, Tozzi & Rosati 
(2003, Paper I). The main differences with respect to Paper I 
are the following: (1) the present dataset contains 60 objects that 
spans homogeneously the redshift range between 0.06 and 1.27, 
whereas in the previous analysis 17 objects were selected, 9 at 
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z < 0.1 and 8 at z > 0.7; (2) several assumptions on, e.g., the 
star mass fraction and the depletion parameter are tested against 
different measurements from observations and numerical simu- 
lations; (3) a more extensive analysis is performed to assess the 
systematics that affect our results on the constraints of the cos- 
mological parameters. 

The outline of our work is the following. In Section 2, we de- 
scribe the cosmological framework that allows us to formulate 
the cosmological dependence of the cluster gas mass fraction. 
In Section 3, the dataset is presented and the physical quantities 
under examination are computed. The method based on using 
the gas mass fraction as a cosmological probe is discussed in 
Section 4. The results on the cosmological parameters are shown 
in Section 5, along with a discussion on their robustness against 
systematic effects. The summary of our findings and the conclu- 
sions are drawn in Section 6. Throughout this work, if not oth- 
erwise stated, we plot and tabulate values, with errors quoted at 
the 68.3 per cent (Icr) level of confidence, that are estimated by 
assuming iJp = TOh^^l km s^^ Mpc^^ and 57,„ = 1 — — 0.3. 

2. The cosmological framework 

We refer to f2b as the baryon matter density, to ili„ as the total 
matter density (i.e. Om = ilb + f^c, where fie is the cold dark 
matter component), to as the dark energy density, that we 
consider both in its static and homogeneous form as a cosmolog- 
ical constant and with an equation of state varying with cosmic 
time as "quintessence" (e.g. Turner & White 1997, Caldwell et 
al. 1998). All these densities are expressed in units of the crit- 
ical density, pc = 3Hq / (SttG), where Hq is Hubble's con- 
stant and G is the gravitational constant. In our computation, 
we neglect (i) the energy associated with the cosmic radiation, 

w 4.16 X 10-5(rcMB/2.726ii:)4, and (ii) any possible 
contributions from light neutrinos, ili,h jq — ^mi,/45.5eV, 
that is expected to be less than 0.01 for a total mass in neutri- 
nos, J2 nil'' lower than 0.62 eV (see, e.g., Hannestad & Raffelt 
2006). Therefore, we adopt the Einstein equation in the form 

+ r^A + f^k — 1, where Slk accounts for the curvature of 
space. From, e.g., Carroll, Press & Turner (1992, cf. eq. 25), we 
can then write the angular diameter distance as 

«ang - (1+^)2 - Hoil+z) |nk|i/2' 

where c?ium is the luminosity distance, S{uj) is sinh(tj), ui, 
sin(w) for fik greater than, equal to and less than 0, respectively. 
In addition, we define 

Eiz) = [an(l + 2)-Vr!k(l + z)' + ^^AA(z)]'^' 

A(.) = exp(3£i±^dz). (2) 

The above equations (i) do not include the evolution with red- 
shift of the radiation component, ilr(l + 2)^, that at z « 1 is 
about 1.4 X 10^"^ and therefore negligible in the overall budget, 
(ii) consider the dependence upon the ratio w between the pres- 
sure and the energy density in the equation of state of the dark 
energy component (Caldwell, Dave & Steinhardt 1998, Wang 
& Steinhardt 1998). Hereafter we consider both a pressure-to- 
density ratio w constant in time that implies A(z) = (1 + z)^^^"^ 
and a simple parameterization of its evolution with redshift (e.g. 
Rapetti et al. 2005, Firmani et al. 2005) 

wiz)^wo + wi— — . (3) 



In particular, the case of a cosmological constant requires w = 

-1. 

3. The dataset 

We consider the sample of 52 galaxy clusters at z > 0.3 pre- 
sented in Balestra et al. (2007; see Table [T]!. The preparation, 
reduction and analysis of this dataset is described below. 

As local measurements, we consider the sample of 8 objects 
at a redshift between 0.06 and 0.23 with a gas temperature higher 
than 4 keV described in Vikhlinin et al. (2006; see Table|2]i. We 
use the best-fit results of the gas temperature and density profiles 
obtained by adopting the same functional form that is applied 
in the present analysis to reproduce the ICM properties of the 
high— z sample. The gas and total mass profiles are then recov- 
ered for both the local and the z > 0.3 sample with the same 
method described at the end of this section. 

The histogram of the redshift distribution of the two samples 
is shown in Fig.[T] 

The only criterion adopted for the selection of these objects 
is the availability of good exposures with Chandra of hot, mas- 
sive, relaxed galaxy clusters over a significant redshift range. No 
selection effect is expected to occur in the application of the gas 
mass fraction method once the clusters are selected to ensure (i) 
the use of the hydrostatic equilibrium equation to recover the to- 
tal mass (round, relaxed objects are required for that), and (ii) a 
negligible contribution from non-gravitational energy in the re- 
gion of interest to allow the use of cross-calibration with numer- 
ical simulations, such as the estimate of the depletion parameter 
(the selection of hot, massive systems dominated energetically 
by gravitational collapse satisfies this condition). 

We have considered only Chandra exposures of the clusters 
listed in Table [T] Data reduction is performed using the CIAO 
3.2 software package with a recent version of the Calibration 
Database (CALDB 3.0.0) including the correction for the de- 
graded effective area of ACIS-I chips due to material accumu- 
lated on the ACIS optical blocking filter at the epoch of the ob- 
servation. We also apply the time-dependent gain correctiorQ 
which is necessary to adjust the "effective gains", which drifts 
with time due to an increasing charge transfer inefficiency (CTI). 
The detailed procedure of data reduction is described in Balestra 
et al. (2007). 

A single thermal component (XSPEC model mekal; 
Arnaud 1996) absorbed by the Galactic column density is fit- 
ted to the spectrum extracted from a circular region with radius 
i?spoc (see Table [U around the cluster chosen to maximize the 
signal-to-noise ratio once the contaminating point sources are 
masked. The only free parameters in the spectral fit are the tem- 
perature, the normalization and the metallicity with respect to the 
abundance table from Grevesse & Sauval (1998) here assumed 
as the reference. These temperature measurements are on aver- 
age 3 per cent higher than the values obtained by assuming the 
standard metallicity table of Anders & Grevesse (1989). 

We assume that galaxy clusters are spherically symmetric 
gravitationally bound systems. We deproject the X-ray surface 
brightness to obtain the electron density profile, no(r), by cor- 
recting the emissivity by the contribution from the outer shells 
moving inwards. We consider the profiles from the deprojection 
of the background-subtracted surface brightness S, with eiTor 
es, up to the radial limit, i?spat, beyond which (S/es) < 2. 
The gas density profiles, normalized to the critical density at the 
cluster's redshift, are shown in Fig. |2] We note that the average 

' http://asc.harvard.edu/ciao/threads/acistimegain/ 
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Table 2. Properties of the local cluster sample from Vikhlinin 
et al. (2006). Column (1): name; column (2): adopted redshift; 
column (3); best-fit spectral temperature; column (4): -R500; 
column (5): gas mass fraction within i?5oo- A cosmology of 
{Ho,n^,nA) = (TOkms^i Mpc^i, 0.3, 0.7) is adopted here. 



Cluster 



keV 



-R500 
kpc 



/gas 

(< -Rsoo) 



A133 

A 1795 

A2029 

A478 

A1413 

A907 

A383 

A2390 



0.057 
0.062 
0.078 
0.088 
0.143 
0.160 
0.188 
0.230 



4.15 ±0.07 
6.10 ±0.06 
8.46 ± 0.09 

7.95 ±0.14 
7.38 ±0.12 

5.96 ±0.08 
4.80 ±0.12 
8.90 ±0.17 



1027 ±38 
1320 ± 47 
1419 ± 32 
1398 ± 62 
1377 ± 49 
1149 ± 32 
983 ± 34 
1489 ± 50 



0.087 
0.108 
0.128 
0.125 
0.112 
0.129 
0.129 
0.147 



± 0.006 
± 0.006 
± 0.007 
±0.011 
± 0.007 
± 0.006 
± 0.007 
± 0.009 



60 obj, 57 with kT>4 keV 



-I— p-i— I— I— p-i— I— I— p-i— I— I— p 



T — I — I — I — I — r- 
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Fig. 1. Redshift distribution of the clusters in the whole sample 
(dashed region) and of the selected hottest objects (shaded re- 
gions). 



profile measured in local (z < 0.3) clusters is definitely more 
peaked than the ones observed at higher redshift, with a cen- 
tral gas density that is a factor of 2 higher than in systems at 
0.3 < z < 0.6 and almost a factor of 5 higher than in objects 
at z > 0.6. This follows recent evidence that cooling cores, at 
relativly higher central density, preferentially form later in time 
(e.g. Santos et al. 2008, Ettori & Brighenti 2008). 

To estimate the gas and total mass profiles, we fit the electron 
density profiles with a functional form adapted from Vikhlinin et 
al. (2005): 



Xi 



0-2X0 
r/ao 
r/ai. 



1 



2\ —1.504+03/2 



4) 



X + )""'"'"t4) 




Fig. 2. Radial gas density profiles in units of the critical density 
at the cluster redshift. Green, black and red lines indicate objects 
at redshift < 0.3, 0.3 — 0.6 and > 0.6, respectively. The bottom 
panel shows the ratio (and relative dispersion) of the mean den- 
sity profiles observed in the redshift range 0.3 — 0.6 (green fine) 
and z > 0.6 (red line) with respect to the local profile. 



The number of free parameters is reduced according to the num- 
ber of datapoints A^dat available, in the following order of prior- 
ity: 05 = 3 (always), ag — 1.2 (when 6 < A^dat < 8, otherwise 
free to vary between and 2.5), 04 — 0.6 (when A^dat — 5), 
03 = (when A'dat < 4). The best-fit parameters that reproduce 
the electron density profiles up to the outer radial limit i?spat 
are quoted in Table [1] All the electron density profiles and the 
best-fit lines are shown in Fig. lA.ll 

The total and gas masses are then evaluated at i?A = ^500 
that describes the sphere within which the cluster overdensity 
with respect to the critical density is A = 500. This value of 
overdensity is maintained fixed in all the cosmologies studied 
and at any redshift. Under the assumption of isothermahty, we 
calculate therefore 



Mtot(<r) = 
-R500 = 

Mgas(<r) = 



firriuG dlogr ' 



3 A/tot (< R. 



•500 J 



47r500/9c 



1/3 



H50 



1.155mu?ic(r) Anr^dr, 



(5) 



where 1.155 is the value associated to a cosmic mix of hydro- 
gen and helium with 0.3 times solar abundance in the remaining 
elements with a relative contribution that follows Grevesse & 
Sauval (1998), /i = 0.600 is the corresponding mean molecular 
weight, TTiu — 1.66 x lO^^"* g is the atomic mass unit and 
is the critical density at redshift z and is equal to 3H'^ / {8ttG) 
with Hz = HqE{z) (see eq. |2]i. The gas mass fraction is then 
/gas(i?A) = Mgas(< i?A)/Mtot(< ^^a)- The distribution of 
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Fig. 3. Ratios between the maximum radial extension in the 
gas density profile, -Rspat, and i?5oo (diamonds) and the extrac- 
tion radius of the circular region used to map the object in the 
spectral analysis, i?spoc, and i?5oo (dots). The dashed and solid 
lines show the average values of these ratios, i?spat/^500 and 
^spccZ-Rsoo respectively, in bins of ^^10 objects. 



^spat / -R500 and of i?spcc 7^500 as a function of redshift is shown 
in Fig. [3] 

All the eiTors are determined at Icr confidence level from 
the 16th and 84th percentile of the distribution of the values ob- 
tained by repeating the calculations 1000 times after a new sur- 
face brightness profile and gas temperature are considered by 
selecting normally-distributed random values according to the 
original measurements and their relative errors. A new total mass 
profile is then estimated and new i?5oo and gas mass measure- 
ments are obtained by following equations |5] 

Our estimates of the gas fraction rely upon the determina- 
tion of the total gravitating mass obtained under the assumption 
of an isothermal gas. To check the effect of the presence of a 
temperature profile on the measured total mass and, hence, gas 
fraction, we assume the expression for the temperature profile 
that Vikhlinin et al. (2006) showed to reproduce well the tem- 
perature gradients in nearby relaxed systems: 



T(r) = 1.23 



(x/0.045) 



1.9 



0.45 



(x/0.045)i' 



1 



[{x/o.ey + 1] 



0.45 ■ 



(6) 



where x = r/R^oQ and the factor 1.23 comes from equations 8 
and 9 in Vikhlinin et al. (2006). Starting from the value of i?5oo 
measured under the assumption of isothermality, we calculate (i) 
the temperature gradient, (ii) the total mass using the hydrostatic 
equilibrium equation (Aftot oc T(r) x r x d log(nc T)/(i log r) 
and (iii) a new estimate of i?5oo- The calculations are repeated 
until i?5oo, that is evaluated from the total mass profiles and 
adopted in equation|6] converges. This happens after 2-6 itera- 
tions. At the newly estimated i?5oo, we measure the gas and total 
masses. We find that introducing a temperature profile decreases 
^500 by 7 per cent, on average, and increases /gas(< -R500) by 



about 16 per cent with respect to the estimates obtained under 
the assumption of isothermality. The changes induced in the de- 
termination of i?5oo and of /gas are shown in Fig.|4l 

As a reference, for (/170, = 1 — ^^a) = (1, 0.3), in the 
overall sample of 52 objects at 2: > 0.3, we measure /gas. 500 — 
(0.106, 0.113, 0.044) (median / mean / dispersion), when a gas 
temperature profile is assumed, and (0.096, 0.099, 0.043), when 
a constant kT is adopted. 

In the following analysis, we make use of the estimates ob- 
tained with a temperature profile to limit the cosmological pa- 
rameters for the case of reference, and leave to the discussion on 
the systematics the case with a constant temperature. 

3. 1. Comparison with previous work 

We compare our results with recent determinations of the gas 
mass fraction obtained from a similar X-ray dataset. 

Allen et al. (2008) presented an extensive cosmological 
study which uses the gas mass fraction in 42 hot (kT > 4 
keV) relaxed clusters. When we consider the 16 objects in com- 
mon with the present work, the mean (dispersion) of the ra- 
tios of the estimated /gas at i?2500 between their measurements 
and our estimates obtained by assuming a temperature profile is 
1.14(0.24). Also within the scatter observed in this small sam- 
ple used for comparison is the ratio between the spectroscopic 
measurements of the gas temperature (TAUcn/Jthiswork — 1.06, 
r.m.s = 0.17). 

We recall here that our gas temperatures are measured in 
spectra obtained to maximize the signal-to-noise ratio, includ- 
ing the contribution from the core emission. Any cool compo- 
nent in the core is expected to bias downwards the emission- 
weighted cluster temperature. We have checked the statistical 
effect of this bias by considering the sample of 90 clusters with 
temperatures > 4 keV for which Maughan et al. (2007) obtained 
from Chandra observations a measure of the temperature both 
in the range (0 — l)i?5oo and in the range (0.15 ~ l)i?5oo- We 
obtain a mean To is-i/Tq-i of 1.01, with a dispersion of 0.14, 
suggesting that, even though on average a cool core does not 
significantly affect the cluster temperature, the scatter in the dis- 
tribution of the measured temperatures can be high. 

Overall, these results suggest that some systematics (e.g., the 
background subtraction, the energy range adopted for spectral 
fitting; see discussion in Appendix of Balestra et al. 2007), of 
the order of 10-20 per cent and comparable to the mean statisti- 
cal error of 10 per cent, affect the estimates of the global temper- 
atures. These uncertainties propagate also to the measurements 
of the gas mass fraction. To assess the role of these eiTors on 
the robustness of our constraints, we investigate in Sect. l5.4l how 
our results depend upon the sample selected according to the 
measured relative error on the gas mass fraction, an error that is 
propagated from the uncertainties on Tgas. 



4. The cluster gas fraction as cosmological 
probe: the method 

We start with a grid of values for the cosmological parameters 
we want to investigate, specifically 57,„ and JIa or w. We com- 
pute then 



x' = 



E 



ifha.r,i/bi — ilb/f^m)^ 



In the above equation, we use the following definitions: 



(7) 
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Fig. 4. Changes in estimates of i?5oo and /gas once the temperature profile in eq.|6]is assumed. Median relative error bars are shown. 



/bar = /gas + /cold, where the gas mass fraction /gas 
is directly measured from our X-ray observation and de- 
pends upon the cosmological parameters through the an- 
gular diameter distance, dang, defined in equation [T] be- 



Mgas/Mtot OC ngasi?^/i? 



oc 



j5/2 



/da 



OC 



ing /gas 

dang(f^m, f^A, wf"^"^ (scc Fig.|5]for the relative dependence 
upon different cosmologies), while the mass fraction in cold 
baryons, /coid — ^^coid/^^tot, is estimated statistically as 
described in the next subsection; 
2. the error on /bar, Cbar, is the sum in quadrature of the uncer- 
tainties on /gas, /star and on the assumed value of Hubble's 
constant (see below) propagated through the following 



dependence: /gas esc iJg 



and /cold (X 



3. the depletion parameter h (with error £{,) represents the frac- 
tion of cosmic baryons that fall in the cluster dark matter 
halo and is estimated from hydrodynamical simulations as 
discussed in section 4.2; 

4. we assume f7b/i?o = 0.0462 ± 0.0012 (eiTor at Icr level) 
from the best-fit results of the joint analysis in Komatsu 
et al. (2008) of (i) the power spectrum of the tempera- 
ture anisotropy measured from the Wilkinson Microwave 
Anisotropy Probe five year data release (WMAP 5-year, 
Dunkley et al. 2008), (ii) a combined set of magnitudes 
of Type la supernovae (Riess et al. 2007, Astier et al. 
2006, Wood-Vasey et al. 2007), (iii) the Baryon Acoustic 
Oscillations measured in a survey of galaxies at z = 0.2 
and z = 0.35 (Percival et al. 2007). 

We note that the recent compilation of the best-fit results 
from the Primordial Nucleosythesis calculations on fib pro- 
vides comparable constraints (see, e.g., Steigman 2006, page 
34: f7b/i?o = 0.0454 ± 0.0045). We consider the WMAP 5- 



year limits in the following analysis. 

We adopt a present-day Hubble's constant of Hq 



72 ± 8 



km s^ (error at \a level) measured in the Hubble Space 
Telescope Key Project by Freedman et al. (2001). This value, 
assumed also as a Gaussian prior in the joint cosmological 
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Fig. 5. Sensitivity of the cluster baryon fraction method to the 
variation of the cosmological parameters. The shaded region in- 
dicates the redshift range considered in the present study. 



analysis from Komatsu et al. (2008), is in agreement with 
their final result of i7o = 70.1±1.3 km s^^, but is marginally 
in conflict with the value of iJo ~ 62.3 ± 1.3(random) ± 
5.0 (systematic) km s^^ Mpc^^ recently determined from 
Type la supernovae calibrated with Cepheid variables in the 
nearby galaxies that hosted them (Sandage et al. 2006). We 



6 



S. Ettori et al.: Cosmological constraints from the Cluster Gas Fraction 



0.030 
0.025 

0.020 

J 0.015 

0.010 

0.005 

0.15 

o 

^0.10 

o 

0.05 



0.00 




<> <> 




0.0 0.2 0.4 0.6 0.8 1.0 
redshift 



1.2 



Fig. 6. Stellar mass fractions (cf. equation[8]l and ratios between 
their values and the estimated gas mass fraction at -R500 as 
function of redshift. Dashed lines indicate the median value: 
/cold — 0.009 and /coid//gas = 0.091. The two diamonds indi- 
cate the values for RDCS-J0910 and RDCS-J1252 with the stel- 
lar masses estimated from the near-infrared luminosity function 
in Strazzullo et al. (2006). A cosmology of {Hq, D,^, Qa) = (70 
km s~^ Mpc^^, 0.3, 0.7) is adopted here. 



make use of these estimates of Hq in section 5, where we 
discuss the robustness of our cosmological constraints. 

For a set of cosmological parameters {ftm, ^a, w}, the gas 
mass, the total mass and the critical density pc^z are evaluated 
by considering all their cosmological dependencies. A new gas 
fraction for each cluster in the selected sample is then estimated 
at the fixed overdensity of A = 500 following equations |5] 



4. 1. The stellar mass fraction 

The stellar mass fraction, /star — ^'^star/Aftot, cannot be com- 
puted individually for each single object for which most of the 
information of, e.g., the luminosity function of many clusters 
at medium and high redshift is lacking. Therefore, a statistical 
approach is generally used by estimating /star as a function of 
other observed quantities for which it has been possible to con- 
firm a robust correlation on a small and local sample of galaxy 
clusters. 

From the measurements of the optical luminosity and Rosat 
X-ray surface brightness in Coma, White et al. (1993) found 
/star = 0.18(±0.05)/gas. Arnaud et al. (1992) quote /star = 
0.17(±0.03)/gas when their values are converted using a mass- 
to-light ratio for stars in the V band of Q{M/L)q (see also 
Voevodkin & Vikhlinin 2004). Similar results have been ob- 
tained from Fukugita, Hogan & Peebles (1998) by estimating 
the global budget of cosmic baryons. In the following analy- 
sis, we will refer to these very similar recipes as /star w = 

0.18(±0.05)/gas. 



Lin, Mohr & Stanford (2003) use the near-infrared light as a 
tracer of the total stellar mass within cluster galaxies. They find 
a good correlation between the stellar mass fraction obtained 
from the K-band luminosity function measured using the Two 
Micron All Sky Survey (2MASS) data and X-ray properties of 
27 nearby systems. Considering that their results apply to all the 
sample that spans a gas temperature between 2 and 9 keV and 
that they rely on total mass estimates obtained from scaling re- 
lations taken from the literature, we use their measurements of 
Afstar and their assumed Tgas to recover an empirical relation 
valid for the clusters with Tgas > 4 keV that we will consider in 
our analysis. We obtain that, for these 10 very massive clusters, 
the ratio Afstar, 500 /^gag is almost independent of the tempera- 
ture itself and is equal to 5.09(±1.73) x lO^^^© (5 keV)-!-^ 
(weighted mean and dispersion after propagation of the errors 
on both Afstar. 500' for which we assume the relative error on 
the estimates of the luminosity, and Tgas). The estimated stellar 
masses are then added to the measured gas masses to evaluate the 
total baryon mass. We will refer to /star,L to indicate this method 
calculating the stellar contribution to the total baryon budget. 

An additional component of the cold baryon budget is the 
intracluster light (ICL) at very low surface brightness, mak- 
ing the total cold baryonic content of galaxy clusters equal to 
/cold — /star + /iCL- Numerical simulations (e.g. Willman et 
al. 2004, Murante et al. 2004, 2007) suggest that the ICL mostly 
originates from the merging processes taking place during the 
assembly of massive cluster galaxies. Consistently with observa- 
tional results for the Virgo cluster (Feldmeier et al. 2004), sim- 
ulations predict /icl ~ 0.2/star in lower mass systems, with an 
increasing fraction up to /icl ~ /star in 10^^ A/© clusters (e.g. 
Lin & Mohr 2004). 

Gonzalez, Zaritsky & Zabludoff (2007) have presented a 
census of the baryons in local systems with total masses in 
the range 6 x lO^^Af© - lO^^A/©. Including the ICL luminos- 
ity well mapped within 300 kpc, but probably underestimated 
in the outskirts if there is any significant contribution at larger 
radii, they found a well defined correlation between the stellar 
mass fraction and the total mass at A = 500 in the form of 
/cold - (0.009 ±0.002)(Af5oo/10i5Afo)-(0-64±o-i3). 

Recently, Lagana et al. (2008) discussed in detail the bary- 
onic content of five massive galaxy clusters, including an ICL 
contribution. They conclude that the stellar-to-gas mass ratio 
within i?5oo anti-correlates slightly with the gas temperature and 
can be expressed through the relation 



/ci 



.Id 



/star + /icl = (0.18 - 0.012Tgas) /g 



(8) 



where Tgas is measured in keV. For clusters with Tc 



> 4 keV, 

that is one of our selection criteria, /coid//gas < 0.12. These val- 
ues lie on the lower end of the stellar mass distribution presented 
in the above-mentioned work, whereas /star,w represents an up- 
per limit. However, considering that the analysis has been per- 
formed self-consistently for the baryonic component for a sam- 
ple of massive galaxy clusters that well fit the properties of the 
objects studied in our work, we adopt this correlation between 
/cold and /gas and propagate the estimated 20% uncertainty on 
the /cold value. 

As a reference, for (/i7o,rim — ^ — ^h) = (1,0.3), we 
obtain for the selected sample of massive objects /coid,500 — 
(0.091/0. 090/0. 027)/gas, 500 (median / mean / dispersion; see 
Fig.|6]l which represents about 10% of the total cluster baryon 
budget (see discussion in Section |53T l. 

We study further the effect of (i) a different recipe to eval- 
uate /star and (ii) the possibility of a larger, up to 20 per cent. 
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Fig. 7. Depletion parameter as a function of redshift as mea- 
sured in cosmological hydrodynamical simulations in Ettori et 
al. (2006). The solid and dashed lines show results for non 
radiative simulations and for radiative simulations, which in- 
clude star formation and galactic winds, respectively. These es- 
timates are averaged over simulations of four massive clusters 
with M200 > IO^^Mq/H at z = 0. The dots and errorbars in- 
dicate the average 6500 and the scatter around it. The red solid 
line indicates the best-fit constant value 6500 estimated at 2; = 
in non-radiative simulations (dotted lines show the dispersion 
around the mean). The green solid line represents 6500,2 of equa- 
tion ( fTOl ) adopted in our analysis. 



contribution of the ICL, when we discuss the robustness of our 
results to some systematic errors in Sect |5.4| 



4.2. The depletion parameter 

The depletion parameter b — /bai/(f^b/f^m) indicates the 
amount of cosmic baryons that are thermalized within the cluster 
potential. When measured over well-representative regions of a 
galaxy cluster, this value tends to be lower than unity because not 
all the accreting shock-heated baryons relax within the dark mat- 
ter halo. However, the depletion parameter can presently only be 
inferred from numerical simulations, where both the input cos- 
mic baryons are known and the amount of them accreting into 
the cluster dark matter halo can be traced. For instance. Eke, 
Navarro & Frenk (1998) find in their smoothed particle hydro- 
dynamics non-radiative simulations that the gas fraction within 
i?vir is, on average, 87 per cent of the cosmic value. On the 
other side, recent work with grid-based numerical codes esti- 
mates that a larger amount of baryons are captured in the clus- 
ter potential (e.g. Kravtsov, Nagai & Vikhlinin 2005). In a set 
of SPH and Eulerian simulations of a single cluster presented 
in the Santa Barbara Comparison Project, Frenk et al. (1999) 
measure at the virial radius h ~ 0.92 ± 0.07. In Ettori et al. 
(2006), we analyze the h parameter in a set of four massive 
(Af2oo > lO^^ft-^^M0 at z = 0) galaxy clusters simulated by 



using the Tree-nSPH code GADGET-2. We consider here the 
distribution of the values of 6(< i?5oo) measured in two sets 
of hydrodynamical simulations, one with gravitational heating 
only, the other including radiative processes, such as cooling 
and star formations with feedback provided from weak winds 
(Fig. 01. We measure h{< i?5oo) = 0.874 ± 0.023 at z = 
and b{< R^qq) = 0.947 ± 0.037 at z = 1 in the first set, 
and 6(< i?5oo) = 0.920 ± 0.026 at z = and b{< R^qq) = 
0.961 ± 0.028 at z = 1, when cooling and feedback are con- 
sidered. If we average over the redshift range 0-1, we obtain 
&(< Rboa) = 0.940±0.023, or, by fitting apolynomial at 1st or- 
der in redshift, 6(< i?5oo) 0.923(±0.006) + 0.032(±0.010)z. 
Considering that the four objects under consideration reached 
^^200 > 10^^ h"^ Mq at z = after they accreted mass by a 
factor 2-6 from z = 1, we adopt two values for the depletion pa- 
rameter, encompassing the results of our simulated dataset (see 
Fig. 13: 

6500 = 0.874 ± 0.023, (9) 

that is related to the most massive systems, subjected only to the 
gravitational heating, in the studied simulation and, thus, par- 
ticularly appropriate for the high-temperature objects that are 
analyzed in the present work over radial regions, beyond their 
cores, where radiative processes (like cooling, star formation, 
stellar feedback) are not expected to influence significantly the 
ICM physiscs; 



"500, z 



0.923(±0.006) + 0.032(±0.010)z, (10) 



that is the quantity measured when all the simulated dataset are 
considered over the redshift range of interest. 



5. Results on the cosmological parameters 

We measure the gas mass fraction at i?5oo in 52 z > 0.3 objects 
with a median relative error tgas of 36 per cent (Icr). In the local 
sample, the relative eiTor is about 6 per cent. 

In the following analysis, we focus on the most massive ob- 
jects present in our sample to work under the ideal hypothesis 
that the physics of the X-ray emitting plasma is only determined 
by the gravitational collapse. To define a sample of massive sys- 
tems whose baryon content is expected to be representative of 
the cosmic budget, we select the 49 objects at z > 0.3 with 
Tgas > 4 keV (see Fig. [1] and Table [1]). We also consider the 
8 galaxy clusters with Tgas > 4 keV in the local sample of 
VikhHnin et al. (2006; see Table|2ll 



5. 1. Constraints on fl,.. and 



We assume hereafter that ilt,h^ 



72 ± 8 km s"^ Mpc" 
(0.18-0.012Tga.)/, 



70 



0.0462 ± 0.0012, Ho = 
6500 = 0.874(±0.023) and /eoid = 
with an associated 20 per cent relative 
error (at la level). From the estimated mean gas mass fraction 
within i?5oo of 0.11 (see Sect. 3), the resulting matter density 
parameter is expected to measure (see eq. |7]) 



bQv 



/gas + /( 



0.32. 



(11) 



cold 



To make extensive use of the information present in the gas 
mass fraction distribution, we minimize eq.|7]for a grid of values 
of {f2m,f2A} that satisfy the condition ft^ + Ha + flk = 1 
(w = —1). We obtain the best-fit results presented in Fig.|8]and 
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Fig. 8. Constraints in the fi„i — Qa plane (w — —1) provided from the sample of 57 objects. (Left) Likelihood contours at 1, 2 and 
3 a {Ax^ = 2.30, 6.17, 11.8, respectively, for 2 degrees of freedom) for the constraints obtained through equation |7] (solid lines) 
and equation [16] by using the gas fraction evolution only (dashed lines). (Middle) Marginalized probability distribution for f2,„. 
(Right) Marginalized probability distribution for Ha- The shaded regions show the constraints from the joint cosmological analysis 
of CMB+Snla+BAO data in Komatsu et al. (2008). 
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Fig. 9. (Left) Constraints in the flm — w plane (fim + f^A = 1) plane provided from the sample of 57 objects. Likelihood contours 
at 1, 2 and 3 cr (Ax^ = 2.30,6.17, 11.8, respectively, for 2 de grees of freedom) for the constraints obtained through equation |7] 
(solid lines) and equation[T6]by using the gas fraction evolution only (dashed lines). (Middle) Marginalized probability distribution 
on w overplotted on the best-fit constraints obtained from the joint cosmological analysis of CMB+Snla+BAO data in Komatsu et 
al. (2008). (Right) Constraints in the wq — wi plane from the definition adopted in equation[3] 



quoted here at the la level of confidence on one single parameter 
(X^ — Xmin — 1) marginalization over the others; 

nA = 0.5910:^6. (12) 

The total is 58.9 for 57 — 2 degrees of freedom (see Fig.fTOli. 

These results change slightly once the recent determina- 
tion of Bubble's constant. Ha = 62.3 ± 1.3(random) ± 
5.0 (systematic) km s^^ Mpc^^ by Sandage et al. (2006), is 
adopted. With a — 66.5, we obtain 



flA = 0.63 



37+0 04 

'-'•'-"-0.04 
+0.46 
-0.48- 



(13) 



More stringent results are obtained when tighter constraints 
on Hubble's constant provided with a joint-fit with the cosmic 



microwave temperature anisotropies are considered. For exam- 
ple, by using the best-fit values from WMAP 5-year data analy- 
sis by Komatsu et al. (2008; Hq = 70.1 ± 1.3 km s^^ Mpc~^, 
f^b = 0.0462 ± 0.0015; labelled "WMAP5" in TableO, we ob- 
tain 



r, O9+0.03 
"■''^-0.02 

1 01 +"■2° 

^•"-■^-0.28' 



(14) 



but with ax ~ 148.9 mainly due to the propagation of the 
relative error on the Hubble's constant value that is smaller than 
in our reference case (2 per cent vs. 1 1 per cent). 
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Fig. 10. Distribution of the measured baryon fractions in the 
sample of 57 objects (49 of which are at z > 0.3) with Tgas > 4 
keV. The error bars at la include all the components listed in 
Section 4. The best-fit results of the cosmological model used as 
the reference are assumed. In the lower panel, the residuals with 
respect to the assumed value of fhai/b x flj^/^h = 1 are shown 
(solid line). The dashed line refers to the distribution obtained 
by fixing all the cosmological parameters to the best-fit results 
and assuming J^a = 0. 



5.2. Constraints on the equation of state oftlie dark 
energy 

We investigate the constraints that the cluster baryon fraction 
can place on the equation of state of the dark energy, w — P/p, 
where P and p are the pressure and density of the dark energy 
component, respectively. We consider both the case w ^constant 
and w — 'w{wo,wi) as in equation[3] We present our reference 
results in Fig. |9l 

Under the assumption of a flat Universe (fim + = 1), the 
best-fit values for our reference model are (x^ = 58.9; la level 
of confidence on one single parameter): 



n 00+0.04 
'-'•'-'^-0.05 



w 



-1.1 



+0.60 
-0.45- 



(15) 



No significant limits are obtained for an evolving w (see 
eq.O, with wq — — I.O^q'y and no constraints at 1 ct level on 
wi in the range (—3, 3). 

Similarly to the case for JIa, changes in the assumed prior 
distributions induce appreciable differences in the central values 
of the best-fit parameters. As shown in Fig.[TTl these differences 
are however always within 2a for w, even in the most extreme 
cases, e.g. not considering the local (z < 0.3) objects or includ- 
ing in the analysis only the hottest (kT > 8 keV) clusters. For 
the measurements of the behaviour is the same as discussed 
in the previous section. 



5.3. Do ttie data require f^A > ? 

Considering our best-fit results presented in Section ISTl we ob- 
tain that JIa > with a significance in the range 1 — 3.6a (see, 
e.g.. Table O. In particular, for our reference case (Fig. [8]), our 
dataset requires a Ax^ of about 1 when JIa is fixed to zero (see 
Fig.fTOl). Therefore, the baryonic mass fraction in galaxy clusters 
alone can only provide marginal evidence for a no null contribu- 
tion from a dark energy component. This contribution is mainly 
due to the assumption that the gas mass fraction does not evolve 
with redshift in the "correct" cosmology. We can write it as 



E 

1=1 



gas, 2 



4 



(16) 



where /gas.i and egas.i are the single gas mass fraction measure- 
ments with the relative errors, /gas and ej are the mean and the 
error on the mean of the /gas,i values. In Fig. [8] the likelihood 
contours are overplotted on the best-fit results. From these plots, 
it is possible to appreciate the different role played by the two 
assumptions made to define the "correct" cosmology: while the 
condition that the gas fraction is constant with look-back time 
marginally influences the value of fini but, at the same time, sig- 
nificantly tilts the contours favouring an ft^ > for il„i ^ 0.35, 
the condition that the cluster baryonic fraction is representative 
of the cosmic value has the greatest statistical weight, collapsing 
the probability around the best-fit value of flm ~ 0.35. 

5.4. Robustness to systematic errors 

In this section, we investigate how our cosmological constraints 
are affected by potential systematic eiTors. To evaluate this, we 
repeat the analysis described in the previous sections and apply 
it to several different subsets defined by changing one of the cri- 
teria adopted for the selection of the reference sample of 57 hot 
galaxy clusters. 

First of all, we do not observe any dependence of /gas on 
both temperature and metallicity in our reference sample: once 
objects with Tgas > 4 keV are selected, the Spearman's rank 
correlation is about 0. 1 implying an underlying correlation sig- 
nificant at less than the la confidence level. 

Moreover, we notice that the estimated gas mass fraction 
within i?2500 is about 75 per cent of the value measured at i?5oo 
(mean and standard deviation after 1000 bootstrap resampling 
= 75 ± 2; median — 73 per cent), confirming its increased value 
outwards. 

As shown in the plots of Fig.[TT|(and summarized in TableO, 
the largest uncertainties are related to the assumed prior on the 
temperature profile. By using the values obtained by assum- 
ing an isothermal ICM (see /gas in Table [T]), and applying the 
method described in Section 4, we obtain (total — 75.5; 

an + ^^A + ^^fe = i,w = -1) 



r!,n = o.4oj:°:°3 

< 0.46, 
and (x^ = 76.3; 17,„ + r^A = 1) 
O - 29+° "^ 



W < 



-0.8. 



(17) 



(18) 



The total shows a significant increase (Ax^ = 17), suggest- 
ing a poorer representation of the distribution of /gas and point- 
ing out the necessity for a robust and reliable determination of 
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Fig. 11. Best-fit results on f7,n and 57a {left) and Vl^ and w {right) for the different sources of systematics discussed in the text. 
The shaded regions show the constraints from WMAP 5-year results (Komatsu et al. 2008). All the errors are at the la level. The 
number of clusters that meets the selection criteria and the total are indicated in the top panel. 



the temperature profiles in galaxy clusters, possibly out to the re- 
gion where the gas mass fraction measurements are more repre- 
sentative of the cluster baryon content (r <; -R2500 ~ 0.3i?2oo)- 

The assumption on the stellar mass fraction does not sig- 
nificantly affect the estimates of the cosmological parameters. 
As discussed in Section 14.11 the stellar mass fraction that we 
adopt here from Lagana et al. (2008) lies on the lower end of 
the distribution available in the literature. By assuming instead 
/coid//gas = 0.18, independently from the cluster mass, as gen- 
erally done in similar work (e.g. Ettori et al. 2003, Allen et al. 
2008), we measure a decrease in ilni by about 9 per cent (see 
case labelled "/star.w")- 

To select our sample of 57 clusters, we have applied only the 
selection criterion of Tgas > 4 keV to consider only the most 
massive systems for which the ICM physics should be mostly 
determined by the process of gravitational collapse. 

Once we add the further criterion that the relative uncer- 
tainty on the gas mass fraction values has to be lower than 0.5 
(figas/Zgas < 0.5), the sample is reduced to 41 objects and the 
good fit obtained (x^ = 34.6) provides marginally looser con- 
straints on the cosmological parameters, but in better agreement 
with WMAP 5 -year limits. On the other hand, when very hot 
systems (Tgas > 8 keV) are considered, the sample is reduced to 
only 22 objects and estimates of 51a consistent with zero are ob- 
tained. Similar constraints are obtained when the local {z < 0.3) 
8 clusters are excluded due to the insensitivity of the geometrical 
part of the method to the dark energy component at high redshift. 

In our selection, we do not consider the morphological as- 
pect of the cluster. To study this effect, we have considered the 
measurements of the centroid shift for the 43 objects in com- 
mon with the sample of 90 clusters analysed by Maughan et al. 
(2007) with Tgas > 4 keV. We consider as relaxed objects the 
18 clusters with a centroid shift lower than the median value of 
1.18 X 10~^ measured over the entire sample of 115 clusters 
(see their Table 2). These 18 objects are part of the first sample. 



labelled ceriQ. A second sample, ceni, includes the 16 objects 
for which an estimate of the centroid shift is not available, for 
a total number of 34 clusters. Through these two samples, we 
obtain very similar results, with a relativly larger (smaller) con- 
tributions from the dark energy (matter) component compared to 
the reference values. 

Overall, we can attribute a systematic relative error of about 
±5 per cent to the most robust cosmological constraint, ilm- It 
is estimated by considering the standard deviations of the val- 
ues plotted in Fig. [TT] and summarized in Table [3] around our 
best-fit reference value of 51m = 0.35. Similarly, we associate 
a systematic uncertainty of (—0.25, +0.30) to our best-fit value 
51a = 0.59 (about 45 per cent) and a relative systematic uncer- 
tainty of (—25, +34) per cent to our best-fit value of w = —1.1. 



5.5. The reverse situation: tlie baryon distribution in X-ray 
galaxy clusters after WMAP 5-year 

In the previous sections, we have looked at the best-fit cosmo- 
logical parameters that better reproduce the observed distribu- 
tion of the cluster baryonic mass fraction. We now reverse the 
the problem, by fixing the cosmological parameters to the val- 
ues suggested from the analysis of the anistropies in the cosmic 
microwave background and investigating the average physical 
properties of the baryons in our X-ray luminous galaxy clusters. 

Given the above measurements of the cluster gas and total 
masses, and the estimated contribution from the stellar com- 
ponent, we can estimate the relative distribution of baryons in 
galaxy clusters assuming that the cosmology is defined in its 
components from best-fit results obtained after the WMAP 5- 
year release (Dunkley et al. 2008), combined in a joint analysis 
with supernovae Type la data and Baryon Acoustic Oscillations 
constraints (see Komatsu et al. 2008). In the present work, we 



S. Ettori et al.: Cosmological constraints from the Cluster Gas Fraction 



11 




G 




0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 12. (Left) Cluster baryonic pie obtained by fixing the cosmological parameters to the best-fit results after the WMAP 5-year 
release (Komatsu et al. 2008) and considering the mean values of the hot (/gas), cold (/cold) and other (/ob) baryon fractions 
normalized to the cosmic value, fo(rib/f^in)wMAP- Ranges within two standard deviations (95.4 % level of confidence) from the 
mean (obtained after 1000 bootstrap resamplings) are indicated within the parentheses. (Right) Dependence of the average /ob value 
upon the cosmological parameters {ilm, ^a)- The dotted line indicates a flat Universe. 



follow the analysis described by Ettori (2003) and update the 
results discussed there. 

We estimate the relative contribution of the cluster baryons 
to the cosmic value (r2b/i^m)wMAP — 0-165 ± 0.009, where 
the best-fit results after the WMAP 5-year release are consid- 
ered (Hq = 70.1 ± 1.3 km s-i Mpc-\ = 1 - J^a = 
0.279 ± 0.013, r^b = 0.0462 ± 0.0015; see Table 1 in Komatsu 
et al. 2008). We correct our baryon fractions by the depletion 
parameter with consistent results both when we use b =constant 
and b — b{z). In the panel on the left of Fig. [12] we show the 
mean values (after 1000 bootstrap resamplings) of the ratios be- 
tween the mass fractions in hot / cold / the remaining baryons 
and the cosmic value equals to 6 x (rib/f^m)wMAP- 

The requirement for the presence of other baryons apart from 
the hot. X-ray emitting component and the cold, stellar phase 
can be explained by either (i) an underestimate by a factor of 
3 of the mass fraction in the form of cold baryons, or (ii) an 
underestimate by about 20 per cent of /gas. 

Recently, McCarthy, Bower & Balogh (2007) obtained sim- 
ilar results by combining /gas profiles collected from the liter- 
ature, and corrected for (i) stellar mass contributions, (ii) vio- 
lations of the hydrostatic equilibrium equation and (iii) the de- 
pletion parameter, with the WMAP 3-year results (Spergel et al. 
2007), concluding that the most likely explanation is that fim 
must lie in the range 0.28 — 0.47, i.e. up to 70 per cent larger 
than the best-fit from the WMAP 5-year analysis. 

We emphasize, however, that our conclusion on a no null 
contribution from other baryons not accounted for in X-ray emit- 
ting plasma and cold stars and intracluster light is drawn mainly 
from fixing the cosmological parameters to the results of the 
analysis after the WMAP 5 -year release that strongly support 
a flat Universe solution. On the other hand, our best-fit results 
seem to suggest either a higher ilm value or, if flm from the 



WMAP 5-year constraints is adopted, a higher (lower) Oa(w) 
value, lowering the requested amount of /ob (see panel on the 
right in Fig.fT2l). 



Moreover, by assuming that the WMAP 5-year data provide 
the "correct" cosmology and, thus, fixing {Hq, = 1 — J^a) 
to those values, we can attempt to limit the contribution of other 
effects that can further bias the estimates of the gas mass frac- 
tion. Among these, we consider the level of the ICM dumpi- 
ness (e.g. Mathiesen et al. 1999) and the presence of non- 
thermal pressure support (e.g. Dolag & Schindler 2000). To 
weigh their contribution, we look for a minimum in flie pa- 
rameter space {a,C}, where a = Pnoh- thermal /^thermal is 
the ratio, assumed constant in time and space, between a pres- 
sure of non-thermal origin and the pressure of the ICM and 
C — {< rigj^g > / < "^gas >^)°-^ is a measure of the dumpi- 
ness of the ICM which is observed in brightness as an average 
value of the square of the gas density < n1^^ >. These two pa- 
rameters affect the gas mass fraction estimates through propaga- 
tion of corrections on the total and gas masses in a degenerative 
way: Mtot.now = (l + a)Mtot,oid, Mgas.new = Mgas,oid/C and, 
thus, /gasaicw = /gas,oid/[C(l + a)]. The net effect is thus to 
further reduce the gas mass fraction value. Adding this to the in- 
dication that the WMAP 5-year cosmology does not account for 
the observed baryon census will further increase the measured 
discrepancy. Indeed, we obtain that, at a 3a level of confidence 
with 2 degrees of freedom, our gas mass fraction measurements 
within i?5oo in a WMAP 5-year concordant Universe constrain 
a < 0.13 and C < 1.01 (the best-fit results require a — and 
C= 1). 
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6. Conclusions 

We study the baryonic content in 52 X-ray luminous galaxy clus- 
ters observed with Chandra in the redshift range 0.3 - 1.273. 
We include in the sample the measurements of 8 objects with 
kTgas > 4 keV lying in the redshift range 0.06-0.23 presented 
in Vikhlinin et al. (2006). We adopt the same functional form to 
estimate the gas density and temperature profiles to recover the 
gas and total mass radial distribution. We estimate statistically 
the average contribution from cold baryons in stars and intra- 
cluster light to the total baryonic budget of each cluster. By us- 
ing the baryonic content determined in this way, we have inves- 
tigated the robustness of the cosmological constraints provided 
from the cluster baryon fraction alone. 

We show that the a-priori knowledge of the cosmic baryon 
density fib and the Hubble constant Hq are essential to place 
significant limits in the fim ~ JIa/w plane. In particular, the 
determination of their central values are required to avoid the 
introduction of systematic differences in the estimate of the cos- 
mological parameters. This systematic bias cannot be coiTected 
by simply enlarging the error bars on the assumed central val- 
ues, but need a definitive input from other sources (for instance, 
CMB measurements, type la supernovae data or primordial nu- 
cleosynthesis calculations; see, e.g., Allen et al. 2008). 

We find that the gas mass fraction measured in our subsam- 
ple of 57 hot (Tgas > 4 keV) galaxy clusters, in combination 
with a Bubble's constant value of Hq = 72 ± 8 km s~^ Mpc^^ 
and rib^o = 0.0462 ± 0.0012, provides a best-fit result of 
= 0.35to°ol and flA = 0.59toit (l^r error; w = -1) 
and fi,„ = 0.32;°;°^ and w = -l.lto j,, when a flat Universe 
is considered. 

We discuss in detail the systematic effect on our best-fit con- 
straints originating from (i) the stellar mass fraction component, 
expected to be of the order of 0.1 — 0.2/gas, (ii) the redshift de- 
pendence of the depletion parameter b, (iii) the isothermality of 
the ICM. The latter effect, in particular, increases the total by 
17 and shifts the best-fit values by about 15 per cent, pointing out 
the necessity for a robust and reliable determination of the tem- 
perature profiles in galaxy clusters, possibly out to the region 
where the gas mass fraction measurements are more representa- 
tive of the cluster baryon content (r ^ i?25oo ~ 0.3i?2oo)- 

By fixing the cosmological parameters to the values from 
the WMAP 5-year data analysis, we limit the contributions ex- 
pected from non-thermal pressure support and ICM dumpiness 
to be lower than about 10 per cent. Otherwise, there is room 
to accommodate baryons not accounted for either in the X-ray 
emitting plasma or in stars in the order of 18 per cent of the to- 
tal cluster baryon budget. However, this value is lowered to zero 
once a no-flat Universe is allowed, as suggested from the cluster 
gas mass fraction distribution alone. 

Bearing in mind that the /gas method to constrain the cos- 
mological parameters is not self-consistent and needs a strong a- 
priori knowledge of both the Hubble constant value and the cos- 
mic baryon density, we conclude that it offers well proven and 
reliable results on (relative statistical error at 1 cr of '^1 1%, 
with a further systematic effect of ±5%), whereas it presents 
weakness in terms of the assumptions on the astrophysics in- 
volved when applied to investigate the dark energy issue. On the 
other hand, this method provides a complementary and indepen- 
dent constraint of the parameter space investigated from planned 
dark energy experiments (e.g. Rapetti et al. 2008; see also Linder 
2007 for a more pessimistic view) and should be supported for 
this application. 
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Appendix A: Fit of thie deprojected gas density 
profiles 

We show the electron density profiles, obtained from the depro- 
jection of the surface brightness, of the objects in our high— z 
sample. We overplot the best-fitted function form in equation |5] 
fitted over the radial range — i?spat ■ Both the best- fit parameters 
and i?spat are quoted in Table[T] 
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Fig. A.l. Plots of the deprojected electron density profiles with the best-fit model in equation |5]overplotted and the residuals ai — 
{di — mi)/ei shown in the bottom panel, where di, rrii and are the values of the data, model and eiTor on the data at a given radius, 
respectively. Together with the = "^O^/^f — value, the R factor represents a measure of the overall residuals and 

is given by i? = ^ ahs{di — TOz)/ S 'ii- P^i" ^^^^ of ^h^ cases, R < 0.2. 
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Table 1. Properties of the z > 0.3 cluster sample described in Balestraet al. (2007). Column (1): name; column (2): adopted redshift; 
column (3): best-fit gas temperature within i?spcc; column (4): best-fit parameters of the gas density profile in equation|5]fitted over 
the radial range — i?spat; column (5): outer radius i?spat where the signal-to-noise ratio is above 2; column (6): radius i?spcc of 
the circular region adopted in the spectral fit; column (7-8): i?5oo and /gas estimated by assuming the temperature profile in eq.|6j 
column (9-10): i?5oo and /gas estimated by assuming a constant temperature equal to fcTgas. A cosmology of {Ho, fim, ft^) — (70 
km s~^ Mpc~^, 0.3, 0.7) is adopted here. 
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Table 3. Best-fit results and la (Ax^ = 1) errors of the cosmological parameters constrained in our analysis under different 
assumptions. The quoted constraints are obtained for each parameter after marginalization. In Fig. [TT] the same constraints are 
plotted and the number of clusters in the sample and the total ^rs indicated. The single samples are discussed in Sect. 15.11 15.21 
and El 
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